Unconventional states of bosons with synthetic spin-orbit coupling 
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Spin-orbit coupling with bosons gives rise to novel properties that are absent in usual bosonic 
systems. Under very general conditions, the conventional ground state wavefunctions of bosons are 
constrained by the "no-node" theorem to be positive-definite. In contrast, the linear-dependence of 
spin-orbit coupling leads to complex-valued condensate wavefunctions beyond this theorem. In this 
article, we review the study of this class of unconventional Bose-Einstein condensations focusing 
on their topological properties. Both the 2D Rashba and 3D a ■ p-type Weyl spin-orbit couplings 
give rise to Landau-level-like quantization of single-particle levels in the harmonic trap. The inter- 
acting condensates develop the half-quantum vortex structure spontaneously breaking time-reversal 
symmetry and exhibit topological spin textures of the skyrmion type. In particular, the 3D Weyl 
coupling generates topological defects in the quaternionic phase space as an SU(2) generalization of 
the usual U(l) vortices. Rotating spin-orbit coupled condensates exhibit rich vortex structures due 
to the interplay between vorticity and spin texture. In the Mott-insulating states in optical lattices, 
quantum magnetism is characterized by the Dzyaloshinskii-Moriya type exchange interactions. 

PACS numbers: 



I. INTRODUCTION 

Spin-orbit (SO) coupling plays an important role in in- 
terdisciplinary areas of physics. In quantum mechanics, 
SO coupling arises from the relativistic effect as a low 
energy approximation to the Dirac equation. Its semi- 
classic picture is the Thomas precession that electron 
spin moment couples to a velocity-dependent effective 
magnetic field generated by the Lorentz transformation 
of the electric field. In atomic physics, SO coupling con- 
stitutes one of the basic elements to the formation of the 
atomic structures. The development in condensed mat- 
ter physics shows that SO coupling is indispensable for 
important phenomena ranging from spintronics 1 , anoma- 
lous Hall effects^, spin Hall effects^— , to topological 
insulators^. In particular, topological insulators have 
become a major research focus of current condensed mat- 
ter physics. 

Most current studies of SO coupling are considered for 
fermionic systems of electrons. On the other hand, the 
ultra-cold atomic systems have opened up a whole new 
opportunity to explore novel states of matter that are 
not easily accessible in usual condensed matter systems. 
In particular, it currently becomes experimentally pos- 
sible to implement various kinds of SO coupled Hamil- 
tonians in ultracold atomic gases for both fermions and 
bosons^—. The high controllability of these systems 
makes them an ideal platform to explore novel SO cou- 
pled physics with bosons. 

An important property of bosons is the "no-node" 
theorem, which states that in the coordinate represen- 
tation the many-body ground state wavefunctions are 
positive-definite-^. This theorem is valid under very gen- 
eral conditions such as the Laplacian type kinetic en- 
ergy, arbitrary single-particle potentials, and coordinate- 



dependent two-body interactions. It applies to most 
of the known ground states of bosons, including Bose- 
Einstein condensations (BECs), Mott insulators, and 
supersolids. Technically, it indicates that the ground 
state wavefunctions of bosons can be reduced to positive- 
definite distributions, and thus imposes a strong con- 
straint on bosonic states. For example, it rules out 
the possibility of time-reversal (TR) symmetry breaking 
ground states in traditional boson systems. Consider- 
able efforts have been made in exploring unconventional 
BECs beyond the "no-node" theoremiS.. One way is the 
meta-stable state of bosons in high orbital bands in op- 
tical lattices^—, because the "no-node" theorem does 
not apply to excited states. Unconventional BECs with 
complex- valued condensate wavefunctions have been ex- 
perimentally realized 2 ^—. 

The kinetic energies of SO coupled systems are no 
longer Laplacian but linearly depend on momentum, 
which invalidates the necessary conditions for the "no- 
node" theorem as pointed out in Ref. [l8j]. This pro- 
vides another way towards unconventional BECs. For 
instance, the Rashba SO coupled BECs were early in- 
vestigated for both isotropio 2 ^ and anisotropic cases22. 
In the isotropic case, Rashba coupling leads to degen- 
erate single-particle ground states along a ring in mo- 
mentum space whose radius k so is proportional to the 
SO coupling strength. Such a momentum scale is ab- 
sent in usual BECs in which bosons are condensed to the 
zero momentum state, and thus bears certain similari- 
ties to Fermi momentum in fermion systems. If interac- 
tion is spin-independent, the condensates are frustrated 
in the free space at the Hartree-Fock level, and quan- 
tum zero point energy selects a spin-spiral state based 
on the "order-from-disorder" mechanism. Imposing the 
trapping potential further quantizes the motion around 
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the SO ring which leads to Landau level type quanti- 
zation of the single-particle levels. Under interactions, 
condensates spontaneously break TR symmetry exhibit- 
ing topologically non-trivial spin textures of the skyrmion 
type^i. All these features are beyond the framework of 
"no-node" theorem. 

Recently, SO coupled systems with ultra-cold bosons 
have aroused a great deal of research interest both ex- 
perimental and theoretical. Experimentally, pioneered 
by Spielman's group 1 ^— , BECs with SO coupling in the 
anisotropic ID limit have been realized by engineering 
atom-laser interactions through Raman processes^—. 
Condensations at finite momenta and exotic spin dy- 
namics have been observed 1 ^— . Various experimental 
schemes have proposed to realize the isotropic Rashba 
SO couplmg22r— . On the side of theory, the Rashba SO 
coupled bosons have been extensively investigated under 
various conditions, including the exotic spin structures 
in the free space, spin textures in harmonic traps, vortex 
structures in rotating traps, and the SO coupled quan- 
tum magnetism in the Mott-insulating state o 26 ' 26 ' 27 ' 33 " — . 
Furthermore, a recent progress shows that a 3D a-p-type 
SO coupling can also been implemented with atom-laser 
interaction s 44 ' 61 . This is a natural symmetric extension 
of Rashba SO coupling to 3D dubbed Weyl SO coupling 
below due to its similarity to the relativistic Hamilto- 
nian of Weyl fermions 62 . The Weyl SO coupled BECs 
have also been theoretically investigate d 44 ' 63 " — . 

In addition to the ultra-cold atom systems, recent 
progress in condensed matter systems has also provided 
an SO coupled boson system of excitons. Excitons are 
composite objects of conduction electrons and valence 
holes, both of their effective masses are small, thus rel- 
ativistic SO coupling exists in the their center of mass 
motion. The effects of SO coupling on exciton condensa- 
tions have been theoretically investigate d 26 ' 66 , including 
the spin texture formations^ 6 -. An important experiment 
progress has been achieved in Butov's group 67 ' 68 , that 
spin textures in a cold exciton gas have been observed in 
GaAs / AlGaAs coupled quantum wells from the photolu- 
minescence measurement. 

In the rest of this article, we review the current theo- 
retical progress of studying SO coupled bosons including 
both the 2D Rashba and 3D Weyl SO couplings. Our 
emphases will be on the non-trivial topological proper- 
ties which are absent in conventional BECs. The single- 
particle spectra will be reviewed in Section [TTJ They ex- 
hibit a similar structure to the Landau level quantization 
in the sense that the dispersion with angular momentum 
is strongly suppressed by SO coupling s 26 ' 36 ' 39 ' 49 ' 61 ' 63 ' 69 . 
However, a crucial difference from the usual magnetic 
Landau levels is that these SO coupling induced Lan- 
dau levels maintain time-reversal symmetry, and thus 
their topology belongs to the Z 2 class^ 6 -. The interplay 
between interactions and topology gives rise to a vari- 
ety of topological non-trivial condensate configurations 
and spin textures, which will be reviewed in Section IlIII 
In particular, the 3D condensates with the Weyl cou- 



pling exhibit topological defects in the quaternionic phase 
space. It is exciting to find an application of the beautiful 
mathematical concept of quaternions. In Section [TV] we 
review the SO coupled BECs in rotating traps, which are 
subject to both the Abelian vector potential due to the 
Coriolis force and the non- Abelian one from SO coupling. 
The combined effects of the vorticity and spin topology 
lead to rich structures 7 -^— , including half quantum vor- 
tex lattices, multi-domain of plane- wave states, and giant 
vortices. Furthermore, in the strongly correlated Mott- 
insulators, SO coupling effects exhibit in the quantum 
magnetism as the Dzyaloshinskii-Moriya type exchange 
interactions 7 -^—, which will be reviewed in Section [V] 
Conclusions and outlooks are presented in Section IVll 

Due to the rapid increasing literatures and the limit of 
space, we will not cover other important topics, such as 
SO coupled fermions 1 ^^^^ 2 -, the SO coupled dipo- 
lar boson s 46 ' 73 . and the proposals for experimental 
implementations^— . 



II. THE SO COUPLED SINGLE-PARTICLE 
SPECTRA AND THE LANDAU LEVEL 
QUANTIZATION 

We begin with the single-particle properties. Consider 
the following Hamiltonian of 2D two-component atomic 
gases with an artificial Rashba SO coupling defined as 



2D P 2 

H o ' = 7jjj + Vt P (r)- X R (<T xPy -<T yPx ), (1) 



where p — —ihV; the pseudospin components f and J, re- 
fer to two different internal atomic components; Xr is the 
Rashba SO coupling strength with the unit of velocity; 
Vt p (r) = ^Muj 2 r 2 is the external trapping potential, and 
lu is the characteristic frequency of the trap. Another SO 
coupled Hamiltonian will be considered is the 3D Weyl 
SO coupling defined as 



H, 



3D,W 




p 

2M 



+ V tp {r) - \ w a -p, 



(2) 



where Xw is the SO coupling strength. 

Even though we will mostly consider bosons for the 
Hamiltonians of Eq. [T]and Eq. [2l they possess a Kramer- 
type TR symmetry T = ia y C satisfying T 2 = —1, where 
C is the complex conjugate operation. At the single 
particle level, there is no difference between bosons and 
fermions. Both Hamiltonians are rotationally invariant 
but break the inversion symmetry. The 2D Hamiltonian 
Eq. [T] still possesses the reflection symmetry with respect 
to any vertical plane passing the center of the trap. For 
the 3D Hamiltonian Eq. [2j no reflection symmetry exists. 

These two typical types of SO interactions have re- 
ceived a lot of attention recently in the community 
of ultra-cold atoms due to their close connection to 
condense matter physics. There have already been 
great experimental efforts on realizing spin-orbit coupling 
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through atom- light interactions ^ 1 2 ' 14 ' 1 5 ' 56 . In fact, sev- 
eral proposals for experimental implementations of Eq. [T] 
and Eq. [2] have appeared in literature o 28 i 31 i 32 ' 44 i 50 . 

In this section, we review the single-particle properties 
of Eqs. [I] and [2] focusing on their topological properties. 
In Sec. Ill A\ their Berry phase structures in momen- 
tum space are presented. When the quadratic harmonic 
trap potential is imposed, the Landau-level type quanti- 
zation on the energy spectra appears with TR symmetry 
as shown in Sec. Ill Bl This Landau level quantization 
provides a clear way to understand novel phases of bosons 
after turning on interactions. In Sec. Ill CI wavefunctions 
of the lowest Landau levels of Eq. [T]and Eq. [2] are explic- 
itly provided. The topology of these Landau level states 
are reviewed through edge and surface spectra in Sec. 

mm 



A. Berry connections in momentum space 

Both Eq. [T] and Eq. [2] possess non-trivial topology 
in momentum space. Let us begin with the 2D Rashba 
Hamiltonian Eq. [1] in the free space, i.e., V tp — 0. Its 
lowest single-particle states in free space is not located 
at the origin of momentum space but around a ring with 
the radius k so — MXr/H. The spectra read 



Now in the free space without the trap, the lowest en- 
ergy single-particle states are located around a sphere 
in momentum space with the radius also denoted as k so 
with the value of k so = MXw/h, and the spectra are 



ti 2 ( \' 1 



(7) 



where the subscripts ± refer to the helicity eigenvalues 
of the operator a ■ k. The corresponding eigenstates are 
solved as 



\i>s 



k+l 



sin -g-e 



(8) 



where cf>k and Ok are the azimuthal and polar angles of 
k in the spherical coordinates. The Berry connection of 
the positive helicity states ipt , is 



A(k) = ~tan-|< 



(9) 



which is the vector potential for a unit magnetic 
monopole located at the origin of momentum space, and 
e t j >k is the azimuthal direction of k. Defining Bi — 
\tijiFji, the corresponding Berry curvature is B(k) = 
2^7 ik, where e k is the radial direction of k. 



(3) 



where ± refer to the helicity eigenvalues of the opera- 
tor a ■ (k x z). The corresponding two-component spin 
wavefunctions of plane- wave states \tp k ±) are solved as 



(4) 



where <j) k is the azimuthal angle of k in the xy-plane. 

For bosons, the lower energy branch states with a fixed 
helicity are important. The Berry connection A(k) of 
positive helicity states ip+(k) is defined as 



A(k) = {ip k+ \iV k \ip. 



k+j 



2k 



e 4>k I 



where e^ k is the unit vector along the azimuthal direc- 
tion. The Berry curvature F.^ is defined as Fij (k) — 
d ki Aj(k) — d kj Ai(k). For a loop winding around the ori- 
gin k = (0, 0), the Berry phase is 



<j> dk ■ A(k) = 7T. 



(6) 



This is because a two-component spinor after rotating 
360° does not return to itself but acquires a minus sign. 
Consequently, Fij (k) is zero everywhere except contribut- 
ing a 7r-flux at the origin of momentum space. 

Next we consider the 3D generalization of the Rashba 
SO coupling of the a ■ p type, i.e., the Weyl coupling. 



B. Landau-level quantization in the harmonic trap 
from SO couplings 

The SO couplings in Eqs. [T] and [2] introduce a SO 
length scale even in the free space defined as l so = l/k so . 
The physical meaning of l so is as follows: the low energy 
states of Eqs. [T]and[2]are not of long- wave length as usual 
but featured with large magnitude of momentum depend- 
ing on the SO coupling strength. l so is the length scale of 
wavepackets that can be formed by using the low energy 
states on the SO ring of Eq. Q] or the SO sphere of Eq. 
[2] On the other hand, in the typical experimental setup 
with ultra-cold quantum gases, a harmonic trap is used to 
confine atoms. The trapping potential V(r) = ^moj 2 r 2 
introduces another length scale as It — \/h/(muj) as the 



(5) typical system size. The trap energy scale is E tp = hu 



It is useful to define a dimensionless parameter a — 
It /ho to describe the relative strength of SO coupling 
with respect to the trapping potential. Physically, a is 
the number of wavepackets which can be packed in the 
trap length. In the limit of large values of a, the trapping 
potential gives rise to Landau level type quantizations in 
both 2D and 3D spin-orbit coupling system o 26 ' 50 i 69 ' 86 . 

The terminology of Landau levels in this section is gen- 
eralized from the usual 2D magnetic case as: topologi- 
cal single-particle level structures labeled by angular mo- 
mentum quantum numbers with flat or nearly flat spec- 
tra. On open boundaries, Landau levels systems develop 
gapless surface or edge modes which are robust against 
disorders. We will see that the low energy states of both 
Eq. [T]and Eq. [5] satisfy this criterion in the case a 3> 1. 
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1. TR invariant Landau levels from, the 2D Rashba SO 
coupling 

Let us briefly recall the usual 2D Landau level arising 
from the magnetic field. In the symmetric gauge with 
A = jBz x r, its Hamiltonian is simply equivalent to a 
2D harmonic oscillator in a rotating frame as 

(p - ^A) 2 p 2 1 

H 2 dll = — C — L - = — + -MwV - ujL z , 10) 

' 2M 2M 2 v ' 

where L z = xp y — yp x and uj = ^j^. Inside each Lan- 
dau level, the spectra are degenerate with respect to the 
magnetic quantum number m. Non-trivial topology of 
Landau levels comes from the fact that m does not take 
all the integer values. For example, in the lowest Landau 
level, m starts from and runs all the positive integer 
number. This chiral feature is a TR symmetry breaking 
effect due to the magnetic field. 

Next let us consider the 2D Rashba SO coupling of 
Eq. Q]in the limit of a > 1. The physics is most clearly 
illustrated in momentum representation. After projected 
into the low energy sector of the positive helicity states, 
the harmonic potential in momentum space becomes a 
Laplacian subjected to the Berry connection as 

M 

V tp {V k ) = -w\K % -A % )\ (11) 

where A is given in Eq. [5]with a 7r-flux at the origin of the 
2D k x -ky plane. In momentum space, the trapping po- 
tential quantizes the motion on the low energy spin-orbit 
ring with radius k so , and is mapped to a planar rotor 
problem. The moment of inertial in momentum space is 
If. = k 2 a M k where M k is the mass in momentum space 
denned as M k = l/(Mu 2 ), and the angular dispersion of 
energy is E agl (j z ) = h 2 j 2 z /2I k = ^jj 2 E tp . Due to the 

7T-flux phase at k = (0,0), j z is quantized to half-integer 
values. On the other hand, the radial component of the 
trapping potential in momentum space is just the kinetic 
energy for the positive helicity states 

H K = hd k io 2 (k-k so ) 2 . (12) 

For states near the low energy spin-orbit ring, the radial 
motion can be approximated as ID harmonic oscillations, 
and the energy gap remains as fku. Combining the radial 
and angular dispersions together, we arrive at 

EnrJ. * {«>' + ^2 + \i l - o?)}E tp , (13) 

where n r is the radial quantum number, and ^(1 — a 2 ) 
is the constant of zero point energy. 

The degeneracy over angular momentum quantum 
numbers is a main feature of Landau level quantization. 
For the Hamiltonian Eq. [TJ although its spectra Eq. [TS] 
are not exactly fiat with respect to j z , they are strongly 



suppressed at a > 1, thus these low energy levels are 
viewed as Landau levels. The radial quantum number n r 
serves as the Landau level index and the gaps between 
Landau levels are roughly E tp . For states in the n r -th 
Landau level with \ j z \ < V^a, their energies remain lower 
than the bottom of the next Landau level, thus they can 
be viewed as gapped bulk states. Actually, the similari- 
ties of these SO coupled states to Landau levels are more 
than just spectra flatness but their non-trivial topology 
will be explained in Sec. Ill Dl 

2. 3D Landau levels from the Weyl SO coupling 

The 3D Landau level systems are not as well-known 
as the 2D case of Eq. [TO] Recently, a large progress 
has been made in generalizing Eq. [TU] to 3D with ex- 
actly flat energy dispersion s 38 ' 93 . In particular, they can 
be constructed with the full 3D rotation symmetry by 
coupling spin-i fermions with the SU(2) gauge poten- 
tial. The Hamiltonian is equivalent to a 3D harmonic 
oscillator plus SO coupling as 

H 3D x L = ^ + \Mu 2 r 2 -uL-ff. (14) 

Excitingly, the lowest Landau level wavefunctions of 
Eq. Q3] possess elegant analytic properties, satisfying the 
Cauchy-Riemann-Fueter condition of quaternionic ana- 
lyticity. Just like that the complex analyticity is essential 
for the construction of fractional quantum Hall Laughlin 
states, the quaternionic analyticity is expected to play an 
important role in high dimensional fractional topological 
states. These 3D Landau level states preserve both TR 
and parity symmetry. The 3D Landau levels have also 
been generalized to the relativistic Dirac particles 9 ^. 

Now let us come back to the Hamiltonian of Eq. [5] 
with the 3D Weyl SO coupling and a trap potential. The 
parallel analysis to the 2D Rashba case applies. Again, 
in the limit of ct > 1, after the projection into the sector 
of the positive helicity states, the trap potential becomes 
V tp (\7 k ) = \M{iS7 k — A k ) 2 and A k takes the form of a 
magnetic monopole one in Eq. [51 The problem is re- 
duced to a spherical rotor problem in momentum space 
on the low energy SO sphere with the radius k so . The 
monopole structure of the Berry connection quantizes the 
total angular momentum j to half-integer values. Simi- 
larly to the 2D Rashba case, the low energy spectra are 
approximated as 

Enr,^ « {n r + + ~(1 - a 2 )}E tp . (15) 

Again the angular dispersion is strongly suppressed by 
SO coupling at a > 1. These spectra exhibit quasi- 
degeneracy over the 3D angular momentum good quan- 
tum numbers of j and j z , and thus can be viewed as a 
3D Landau level quantization with TR symmetry. The 
length scale of these Landau level states is also the SO 
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length l so . Topological properties of these Landau level 
states will be studied in Sec. Ill Dl 



C. Lowest Landau level wavefunctions and parent 
Hamiltonians 

The Landau level energy spectra of Eq. [13] in 2D and 
Eq. [15] in 3D are not exactly flat but with weak disper- 
sions over angular momentum quantum numbers. Never- 
theless, parent Hamiltonians based on slight modification 
on Eqs. Q] and [2] can be constructed. Their lowest Lan- 
dau level spectra are exactly flat and their wavefunctions 
can be solved analytically as shown in Eqs. [17] and [19] 
below. These wavefunctions maintain TR symmetry but 
break parity. In the limit of a 3> 1 and for Landau level 
states with angular momenta \j z | < a in 2D or j < a in 
3D, the lowest Landau level wavefunctions of Eqs. [T]and 
[2] are well approximated by these expressions. 

For the 2D case, the parent Hamiltonian is just 

rr2D,P tt2D,R j- / 1G x 

H = H Q -loL z <j z , (16) 

where L z = xp y — yp x and the coefficient u is the same 
as the trap frequency. As shown in Ref. 61], its lowest 
Landau level wavefunctions are solved as 

where (f> is the azimuthal angle; j z — m+ 5; J m is the ru- 
th order Bessel function. The lowest Landau level energy 
is exactly flat as E LLL = (1 - \)fiu. 

In the case of a > 1 and for small values of \j z \ ps 
to < a, the decay of the wavefunctions Eq. [T7] is con- 
trolled by the Bessel functions rather than the Gaussian 
factor. Their classic orbit radiuses scale as p c j t « ml so . 
Since L z linearly depends on p, the effect of the L z a z 
term compared to that of the Rashba one is at the or- 
der of p c j z uj/Xii ~ rn/a 2 -C 1. Thus Eq. [16] is simply 
reduced to Eq. [T] whose lowest Landau level wavefunc- 
tions are well approximated by Eq. 1171 In this case, the 
length scale of Landau level states is determined by the 
SO length l so instead of the trap length l T . The reason is 
that these Landau levels are composed from plane-wave 
states with a fixed helicity on the low energy Rashba ring. 
The confining trap further opens the gap at the order of 
Hui between SO coupled Landau levels. 

On the contrary, in the opposite limit, i.e., \j z \ rs 
to 3> a 2 , the L z a z term dominates and the Rashba 
term can be neglected. In this case, Eq. [16] is reduced 
into p 2 /2M + ^Mu> 2 r 2 — u)L z a z with a z conserved. In 
each spin eigen-sector, it is just the usual Landau level 
Hamiltonian in the symmetric gauge with opposite chi- 
ralitics for spin up and down, respectively. Nevertheless, 
at to ^> a 2 , the approximation of the projection into the 
Rashba ring for Eq. [T] is not valid, and the eigenstates 
are no longer Landau levels. For the intermediate values 



of a < \j z \ < a 2 , the physics is a crossover between the 
above two limits. 

Following the same logic, the 3D parent Hamiltonian 
with exactly flat SO coupled Landau levels is 

rr3D,P tt3D,W f -. ,-, s 

H ' — H — ujL ■ a, (18) 

where L = r x p is the 3D orbital angular momentum, 
and the coefficient of the L-ff term is the same as the trap 
frequency. Again, as shown in Ref. [61j . the lowest Lan- 
dau level wavefunctions of Eq. fTS] are solved analytically 
as 

TpaBjjS?) = e~^{ji(k S0 r)Y +Ml! (n r ) + ij l+1 (k 30 r) 

x y 1... (o, ■}. (19) 

where ji is the Z-th order spherical Bessel function; 
Y±,j,l,je >s are the SO coupled spherical harmonics with 
total angular momentum quantum numbers j = I ± g 
and j z , which are composed of the spherical harmon- 
ics Yim and spin-i spinors. These lowest Landau level 
states are degenerate over all the values of (jj z ) with 
e lll = (| _ 

Following the same reasoning as in the 2D case, in the 
limit of a > 1, we can divide the lowest Landau level 
states of Eq. Q13 into three regimes as j < a, j >• a 2 , 
and a < j < a 2 , respectively. At j < a, the classic orbit 
radius scales as t c j/It ~ ^, and thus a ■ L comparing 
with a ■ p is a perturbation at the order of j/a 2 -C 1. In 
this regime, the lowest Landau level wavefunctions of Eq. 
[2] are well approximated by Eq. 1191 On the contrary, in 
the regime of j 3> a 2 , <? • L dominates over a ■ p, thus the 
eigenstates of Eqs. [T51 and [21 are qualitatively different. 
In this case, Eq. [18] is reduced to the 3D Landau level 
Hamiltonian Eq. [T4l 

D. The ^-stability of helical edge and surface 
states 

Non-trivial topology of the 2D Landau level manifests 
from the appearance of robust gapless edge states. The 
classic radius r c of each Landau level state expands as 
to increases. For example, in the lowest Landau level, 

r c = \fr~nlB where Ib = \J 1S ^ ne cyclotron radius. 
With an open boundary, as to becomes large enough, 
states are pushed to the boundary^. Unlike the flat bulk 
spectra, the edge spectra are dispersive, always increasing 
with to, and thus are chiral and robust against external 
perturbations. Each Landau level contributes one branch 
of chiral edge modes. If the system is filled with fermions, 
when chemical potential p lies in the gap between Landau 
levels, the chiral edge states give rise to the quantized 
charge transport. 

For Landau levels of the Rashba SO coupling in Eq. 
[T] a marked difference is that these states are TR invari- 
ant. The angular momentum j z in Eq. 1131 takes all the 
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half-integer values as j z = ±5, ±|, ±(m + i), and 
thus these states are helical instead of chiral. Since the 
system described by Eq. Q] does not possess translation 
symmetry, the usual method of calculating topological 
index based on lattice Bloch wave structures in Brillouin 
zones does not work. 

Nevertheless, the non-trivial topology should exhibit 
on the robustness of edge states. The trap length It can 
be used as the sample size by imposing an open boundary 
condition at r — It- States with \j z \ < a are bulk states 
localized within the region of r < It- States with \j z \ ~ a 
are pushed to the boundary, whose spectra disperse to 
high energy rapidly. For a given energy E lying between 
Landau level gaps, each Landau level with bulk energy 
blow E contributes to a pair of degenerate edge modes 
ijj±j z due to TR symmetry. Nevertheless, these two edge 
modes are Kramer doublets under the TR transforma- 
tion satisfying T 2 = — 1. The celebrated Kane-Mele Z2 
argument for translational invariant systems can be gen- 
eralized to these rotation invariant systems by replacing 
linear momentum with the angular momentum. If a given 
energy cuts odd numbers of helical edge modes, then any 
TR invariant perturbation cannot mix these modes to 
open a gap. Consequently, the topological nature of such 
a system is characterized by the Z2 index. The effective 
edge Hamiltonian can be derived as in the plane-wave 
basis if the edge is considered locally flat 

Hedge = v ( sin 7]{p x fi) ■ a z + cos T][(p x h) ■ z)] (t? ■ n)J , 

(20) 

where n is the local normal direction on the circular edge 
in the 2D plane; v is the linearized velocity of the edge 
modes; rj is a parameter angle depending on the details 
of the systems. Both terms are allowed by rotation sym- 
metry, TR symmetry, and the vertical mirror symmetry 
in this system. Each edge channel is a branch of helical 
one-dimensional Dirac fermion modes. 

Parallel analysis can be applied to the surface state 
of the 3D Hamiltonian Eq. [2] Again, due to the TR 
symmetry, the surface state is helical instead of chiral. 
Its topological class also belongs to Z2. If the surface is 
sufficiently large, and thus can be locally taken as flat, 
we can construct the surface Dirac Hamiltonian by using 
plane-wave basis based on the symmetry analysis as 

H s f c — w| sin?7(p x <?) • n + cosi][p ■ a — (p ■ n)(a ■ n)] j, 

(21) 

where h is the local normal direction to the 2D surface. 
Due to the local 50(2) rotational symmetry around n, 
the in-plane momentum p x and p y cannot couple to a z , 
thus spin has to lie in the xy-plane. The first Rashba 
term favors spin and momentum with a relative angle of 
90°, while the second favors spin aligning with momen- 
tum. For a general value of 77 depending on non-universal 
surface properties, Eq. 1211 determines a relative rotation 
between spin and momentum orientation at the angle of 
77. 



III. TOPOLOGICAL SPIN TEXTURES AND 
THE QUATERNIONIC PHASE DEFECTS IN A 
HARMONIC TRAP 

In this section, we review the unconventional BECs 
with interactions and SO couplings, including both 
Rashba and the 3D Weyl types, in the harmonic trap. 
The 2D Rashba case is presented in Sec. IIII Al The lin- 
ear dependence on momentum of SO coupling invalidates 
the proof of "no-node" theorem. Consequently, a general 
feature of SO coupled BECs is the complex-valued con- 
densate wavefunctions and the spontaneous TR symme- 
try breaking. For the Rashba case, the skyrmion type 
spin textures and half-quantum vortex were predicted 
in the harmonic trap™ Furthermore, due to the Lan- 
dau level structures of single-particle states, rich pat- 
terns of spin textures have been extensively investigated 
in literature a 26 ! 39 ' 49 . A nice introduction of topological 
defects in the ultra-cold atom context can be found in 
Ref. (96[. Recently, exact diagonalization work has fur- 
ther showed evidence of strongly correlated phase with 
non-trivial topology 9 -^. 

Even more interesting physics shows in 3D Weyl SO 
coupling, which will be reviewed in Sec. IIII CI The non- 
trivial topology of the condensate wavefunction is most 
clearly expressed in the quaternionic representatio n 64 ' 94 . 
Quaternions are a natural extension of complex numbers 
as the first discovered non-commutative division algebra, 
which has been widely applied in quantum physics 9 ^^^. 
The condensation wavefunctions exhibit defects in the 
quaternionic phase space as the 3D skyrmions, and the 
corresponding spin density distributions are character- 
ized by non-zero Hopf invariants. 



A. Half-quantum vortices and spin texture 
skyrmions with Rashba SO coupling 

Let us consider a 3D two-component boson system 
with contact spin-independent interactions and with 
Rashba SO coupling in the xy-plane. Since the Rashba 
SO coupling is 2D, interesting spin textures only dis- 
tribute in the xy-plane. For simplicity, the condensate 
is set uniform along the z-direction, then the problem is 
reduced to a 2D Gross-Pitaevskii (GP) equation as 

r h 2 V 2 

j - ~^Tf +iSA fl (V x o'y,a/9 - V v <r x ,ap) +gn(r,4>) 
+ hlLU 2 r 2 }Mr,<P)=E^ a (r,^), (22) 

where ij) a 'B with the subscript a =j, ! are two-component 
condensate wavefunctions; n(r, <p) is the particle density; 
g describes the s-wave scattering interaction. The inter- 
action energy scale is defined as E int — gN /(nl T L z ), 
where L z is the system size along the z-axis, and the 
dimensionless interaction parameter is defined as j3 = 
E int /{hu T )- 
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FIG. 1: (A) The radial density distribution of spin up and 
down components, and the total density distribution in the 
unit of iV = f d 3 r{\ip t (r)\ 2 + \i{>i{r)\ 2 } at a = 2 and /3 = 5. 
(B) The skyrmion type spin texture configuration is plotted 
in the rcz-plane. From Refs. [26l . 

We start with weak SO coupling, a ~ 1, and with 
weak interactions. In this case, the energy of single- 
particle ground state with j z — ±i is well separated 
from other states. If interactions are not strong enough 
to mix the ground level with other levels, the condensate 
wavefunction remains the same symmetry structure car- 
rying j z = |, or, — i, thus bosons condense into one of 
the TR doublets, 

H (r, 0) = ( g {$« ) , (r, *) = ( ) , 

(23) 

where f(r) and g(r) are real radial functions. In the 

non-interacting limit, f(r) ~ Jo(fco r ) e u t and g(r) ~ 

Ji{kor)e il T as shown in Eq. 1171 Repulsive interac- 
tions expand the spatial distributions of f(r) and g(r), 
but the qualitative picture remains. Therefore, one spin 
component stays in the s-state and the other in the p- 
state. This is a half-quantum vortex configuration which 
spontaneously breaks TR symmetry^. 

One possibility is that the condensate wavefunction 
may take linear superpositions of the Kramer doublet 
in Eq. [23] The superposition principle usually does not 
apply due to the non-linearity of the GP equation. Nev- 
ertheless, if the interaction of the GP equation is spin- 
independent, all the linear superpositions of the Kramer 
doublet in Eq. [23] are indeed degenerate. This is an ac- 
cidental degeneracy at the mean-field level which is not 
protected. Quantum fluctuations remove this degener- 
acy as shown in the exact diagonalization calculation in 



Ref. HH and select either one of i/)±i. In other words, 
quantum fluctuations can induce a spin-dependent inter- 
action beyond the mean-field levels. Certainly, we can 
also prepare the initial state with the average j z per par- 
ticle ±i, say, by cooling down from the fully polarized 
spin up or down state, then ip±i will be reached. On the 
other hand, if an additional spin-dependent interaction 
is introduced, 

H' m t = sfj d 3 r(n t (r) - n;(r)) 2 , (24) 

then even the mean-field level degeneracy is removed. In 
this shown in Ref. [5ll |. the condensate wave- 

functions of ip±i will also be selected. 

The spin distribution of a condensate wavefunction is 
expressed as 

S(r,cf>) = r a (r,<P)^Mr^), (25) 

which is known as the 1st Hopf map. Without loss of 
generality, the condensate of ipi is considered, and its 

S(r) is expressed as 

S x (r,4>) — psm2^f(r) cos tf>, S y (r, 4>) = p sin 2^(r) sin tf>, 
S z {r,<t>) = pcos2 7 (r), (26) 

where p(r) = \/\f(r)\ 2 + \g(r)\ 2 , and the parameter an- 
gle j(r) is defined through 

cos 7 (r) = 44> sin 7 (r) = 44- ( 27 ) 
p{r) p(r) 

Since the Fourier components of f(r) and g{r) are lo- 
cated around the Rashba ring in momentum space, they 
oscillate along the radial direction with an approximated 
pitch value of ko as shown in Fig. [TJ(A). Because f(r) and 
g(r) are of the s and p-partial waves, respectively, they 
are with a relative phase shift of J. At r = 0, f(r) is at 
maximum and g(r) is 0. As r increases, roughly speak- 
ing, the zero points of f(r) corresponds to the extrema 
of g(r) and vice versa, thus 7 (r) spirals as r increases. 
At the n-th zero of g{r) denoted r n , 7 (r„) = mr (n > 
and we define vq = 0). 

Consequently, S spirals in the zx-plane along the x- 
axis as shown in Fig. [T]B. The entire distribution of S 
can be obtained through a rotation around the z-axis. 
This is a skyrmion configuration which is a non-singular 
topological defect mapping from the real space R 2 to the 
spin orientation space of the S 2 sphere. If the coordinate 
space is a closed manifold S 2 , this mapping is character- 
ized by the integer valued Pontryagin index iT2(S 2 ) = Z, 
or, the winding number. However, the coordinate space 
is the open R 2 , and p(r) decays exponentially at large 
distance r 3> It, thus the rigorously speaking the cov- 
ering number is not well-defined. Nevertheless, in each 
concentric circle r n < r < r Jl+1 , 7 (r) varies from nix to 
(n + 1)tt, which contributes to the winding number by 1. 
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FIG. 2: The phase boundary of f3 c v.s a between (I) the 
skyrmion condensates with j z = ±i and (II) rotational sym- 
metry breaking condensates. From Ref. 12a . 



main driving force towards BEC. For spin- independent 
interactions, the condensate wavefunctions exhibit de- 
generacy at the Hartree-Fock level regardless of the rel- 
ative weight between these two plane-wave components. 
This is a phenomenon of "frustration" . Quantum zero- 
point energy from the Bogoliubov quasi-particle spectra 
selects an equal weight supposition through the "order- 
from-disorder" mechanism. Such a condensate exhibits 
spin-spiral configuration. 

Various literatures have also studied the case of spin- 
dependent interactions in which the Hartree-Fock theory 
is already enough to select either the spin-spiral state, 
or, a ferromagnetic condensate with a single plane- wave 
componen t 33 ' 34 . 



If we use the trap length scale It as the system size, the 
winding number is roughly at the order of a. 

The radial oscillation of the spin density is in anal- 
ogy to the Friedel oscillations in Fermi systems. Around 
an impurity in electronic systems, the screening charge 
distribution exhibits the radial oscillation on top of the 
enveloping exponential decay. The oscillation pitch is 2k f 
reflecting the discontinuity of the Fermi distribution on 
the spherical Fermi surface. Different from the usual bo- 
son systems, the SO coupled ones have a low energy ring 
structure in momentum space in analogous to the Fermi 
surface, thus in real space spin density also oscillates in 
the presence of spatial inhomogeneity. 

In the regime of intermediate SO coupling strength, 
the level spacing between single particle states within 
the same Landau level is suppressed as shown in Eq. 
IT31 In the case that interactions are strong enough to 
mix energy levels with different angular momenta in the 
same lowest Landau level but not among different Lan- 
dau levels, condensates do not keep rotation symmetry 
any more. The calculated phase boundary of interac- 
tion strength /3 v.s. SO coupling strength a is plotted 
in Fig. [2j In this regime, the distributions are no long 
concentric but split into multi-centers and finally form 
a triangular skyrmion lattice structure as calculated in 
Refs. [49], [39]. This 2D skyrmion lattice structure is a 
characteristic feature brought by SO coupling. 



B. Plane- wave type condensations with Rashba SO 
coupling 

If interactions are strong enough to mix states in dif- 
ferent Landau levels, then the influence of the confin- 
ing trap is negligible. The condensate configurations in 
the free space was calculated beyond the mean-field GP 
equation level in Ref. [2d| . Bosons select the superpo- 
sition of a pair of states with opposite momenta k and 
— k on the low energy Rashba ring to condense. These 
spin eigenstates of these two states are orthogonal, thus 
the condensate can avoid the positive exchange interac- 
tions. As is well known, avoiding exchange energy is the 



C. Quaternionic phase defects of the 3D Weyl SO 
coupling 

Next we review the condensates with the 3D Weyl SO 
coupling in the harmonic trap. The corresponding GP 
equation is very similar to Eq. of the Rashba case. 
Only slight modifications are needed by replacing the 
spatial dimension 2 with 3, and by replacing the Rashba 
term with — ifiA^V • a. Amazingly, in this case conden- 
sate wavefunctions exhibit topological structures in the 
quaternionic representation^ 3 -. 

1. The quaternionic representation 

Just like a pair of real numbers form a complex num- 
ber, the two-component spinor tp = (ipf,ipi) T is mapped 
to a single quaternion following the rule 

£ = + & j + &k, (28) 

where £ = Re0t,£i = Inn/>4.,£ 2 = -Re^j,,£ 3 = Im^t- 
k are the imaginary units satisfying i 2 = j 2 = k 2 — 
ijk = — 1, and the anti-commutation relation ij = —ji = 
k. Quaternion can also be expressed in the exponential 
form as 

£= l£|e" 7 = |e|(cos 7 + W sin 7 ), (29) 

where |£| = + and \£\ 2 = & + £ + & w is the 

unit imaginary unit defined as w = + &j + 

which satisfies lu 2 = — 1; the argument angle 7 is defined 

as COS7 = f /|£| and sin 7 = 

Similarly to the complex phase e 1 ^ which spans a unit 
circle, the quaternionic phases e" 7 span a unit three di- 
mensional sphere S 3 . The spin orientations lie in the S 2 
Bloch sphere. For a quaternionic wavefunction, its cor- 
responding spin distribution is defined through the 1st 
Hopf map defined in Eq. [25] as a mapping S 3 — > S 2 . 
Due to the homotopy group o 101 ' 102 ^(S* 3 ) = Z and 
ir 3 (S 2 ) = Z, both quaternionic condensate wavefunctions 
and spin distributions can be non-trivial. The winding 
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FIG. 3: The distribution of S(r) in a) the zz-plane and in the 
horizontal planes with b) z = and c) z/It = \ with a = 1.5, 
c = 1, and — 30. The color scale shows the magnitude of 
out-plane component S y in a) and S z in b) and c). The 3D 
distribution of S(r) is topologically non-trivial with a non- 
zero Hopf invariant. The length unit in all the figures is It- 
From Ref. H 

number of S 3 — > S 3 is the 3D skyrmion number, and that 
of the S 3 — > S 2 is the Hopf invariant, both are integer- 
valued. 

Let us apply the above analysis to the lowest single- 
particle level with j z = j = | 

^ = x(r ) O) = /(r)F +ii!0 i(n) + i 5 (r)y_ i i 1) i(O), 

(30) 

where Y +M (&) = (1,0) T and Y.^ M (n) = 
(cos9, sm9e l ^) T . As shown in Eq. [T9l in the non- 

interacting limit, f(r) jo(kar)e u t and g(r) w 

ji(kor)e il T . The corresponding quaternionic expres- 
sion is 

^ W M)=p(r)e^« (31) 

where p(r) and "f(r) are defined in the same way as the 
2D case in Eq. [26j the imaginary unit, 

u(Ct) = sin 9 cos 4> i + sin 9 sin 4> j + cos 9 k, (32) 

is along the direction of ft. 



2. Quaternionic phase defects and spin textures with 
non-zero Hopf invariants 

The analysis on the topology of the Weyl condensates 
can be performed in parallel to the above 2D Rashba case. 
Again in the case of weak SO coupling, interactions only 
expand the spatial distribution of f(r) and g(r) in Eq. 1301 




FIG. 4: The distribution of S(f) for parameter values a = 
15, p = 0.8, and c = 1 in horizontal cross-sections with (a) 
z/It = 0, (b) z/It = 0.2, respectively. The spin texture 
lattice exhibits a 6cc configuration as sketched in (c). The 
color scale shows the value of S z in (a) and (b). From Ref. 

m 

from their non-interacting forms. The radial wavefunc- 
tion f(r) and g(r) follow the same oscillating pattern as 
those in the Rashba case, thus so does the parameter an- 
gle j(r) which starts from at the origin and reaches 
mr at the n-th zero of f(r n ) = 0. For the quaternionic 
phase e"( sl - l ' 7 ( r " 1 , its imaginary unit uj(fl) is of one-to-one 
correspondence to every direction in 3D, thus it exhibits 
a non-trivial mapping from the 3D coordinate space R 3 
to S 3 , which is known as a 3D skyrmion configuration. 

For a closed 3-manifold, the Pontryagin index of wind- 
ing number is tt^^S 3 ) — Z, i.e., integer. Again here the 
real space is open. In each concentric spherical shell 
with r n < r < r„ +1 whose thickness is at the scale of 
ho, li r ) spirals from nir to (n + 1)tt, thus this shell 
contributes 1 to the winding number from real space 
to the quaternionic phase manifold. If the system size 
is truncated at the trap length It, again the winding 
number is approximately a. In comparison, in the 2D 
Rashba case reviewed in Sec. IIII A) S(r) exhibits the 2D 
skyrmion configuratio n 26 ' 3 ^ 49 , but condensation wave- 
functions have no well-defined topology due to the fact 
that tt 2 (S 3 ) = 0. 

A comparison can be made with the U(l) vortex in the 
single-component BEC. In 2D, it is a topological defect 
with a singular core. Moving around the circle enclos- 
ing the core, the phase winds from to 2tt, and thus 
the winding number is 1 . The above 3D skyrmion phase 
defect is a natural generalization to the two-component 
case whose phase space is S 3 in the quaternionic repre- 
sentation and is isomorphic to the SU(2) group manifold. 
These 3D skyrmions are non-singular defects similar to a 
ID ring of rotating BEC which carries a non-zero phase 
winding number but the vortex core lies outside the ring. 
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The non-trivial topology of the condensate wavefunc- 
tions leads to a topologically non-trivial distribution of 
spin density S(r). The 1st Hopf map defined in Eq. [23] 
becomes very elegant in the quaternionic representation 
as 



Min Max -it 



1 



S x i + S y j + S z k = 



(33) 



where £ = £o — ^2j — £3^ is the quaternionic conjugate 
of £. For the condensate wavefunction Eq. [3lJ S(r) is 
calculated as 



S x (f) 
S y (r) 

S,(f) 



COS ( 

sine 



g(r) sin 
/ 2 (r)+ g 2 (r)cos2#, 



- smc 
cos<b 



g(r) cos f 
f(r) 



(34) 



which exhibits a perfect axial symmetry around the z- 
axis. S(r) is plotted in Fig. [3] at different cross sections. 
In the xy-plane, it exhibits a 2D skyrmion pattern, whose 
in-plane components are along the tangential direction. 
As the horizontal cross-section shifted along the z-axis, 
S(f) remains the 2D skyrmion-like, but its in-plane com- 
ponents are twisted around the z-axis. According to the 
sign of the interception Zq , the twist is clockwise or anti- 
clockwise, respectively. This 3D distribution pattern of 
S(f) is characterized by an integer valued Hopf invariant 
characterized by ^(S* 2 ) = Z. 

As SO coupling strength increases, condensates break 
rotational symmetry by mixing different states with dif- 
ferent values of j in the lowest Landau level. In par- 
ticular, in the case of a > 1, a 3D lattice structure of 
topological defects is formed as shown in Fig. H) Each 
site is a single skyrmion for the condensate wavefunc- 
tion whose spin distribution is a texture with ap- 
proximately a unit Hopf invariant. Spin textures form a 
square lattice in the ccy-plane with the lattice constant 
d ~ 2irl so = 2-kItIol. As for the horizontal cross-section 
planes z ~ id/2, spin textures also form square lattices 
whose sites sit at the plaquette centers of the former lat- 
tice. This is a body-centered cubic lattice configuration, 
which is a generalization of the 2D skyrmion lattice con- 
figuration in Refs. [49( into 3D. Again, if interaction is 
very strong to mix states in different Landau levels, the 
condensate will become plane- wave- like or superpositions 
of SO coupled plane- waves^. 



IV. VORTEX CONFIGURATIONS OF SO 
COUPLED BECS IN A ROTATING TRAP 

Next we review the vortex configurations of SO cou- 
pled unconventional BECs in rotating traps. From a 
more general framework, the above-considered SO cou- 
pling can be viewed as particles subject to non-Abelian 
gauge fields. On the other hand, the Coriolis force from 
rotation behaves as an effective Abelian vector potential. 
Therefore, in a rotating trap, atom-laser coupling pro- 
vides an elegant way to study the effects of these two 




FIG. 5: From left to right: the density and phase profiles of 
ipt(r) and i>i(r) with parameter values of a — 0.5, /3 = 10, 
p = 0.97. From (a)-(g), 7 is taken as 0.5, 0.25, 0.1, 0.0, -0.1, 
—0.3, and —0.5, respectively. At small values of |^| in [c, d, 
and e], a skyrmion lattice is formed near the trap center. By 
increasing the magnitude of ["y| [ (b) and (/)], the skyrmion 
lattice evolves to the normal vortex lattice. For the large 
value of I — y I = 0.5 [(a) and (g)], the condensates show a lattice 
configuration around a ring. The black circle with an arrow 
indicates the direction of the circulation around the vortex 
core. The unit of length for the figures is It- From Ref. [74| . 



different effective gauge fields. We only consider the ro- 
tating systems with the Rashba SO coupling. 



A. Hamiltonians of SO coupled bosons in a 
rotating trap 

Ultracold atoms in a rotating trap share similar physics 
of electrons subject to magnetic fields due to the similar- 
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ity between Lorentz and Coriolis forces. Depending on 
the experimental implementations of rotation, Hamilto- 
nians can be of different type o 70 ' 71 ' 74 . As pointed out in 
Ref. (7pj . because the current experiment setup breaks 
rotation symmetry, rotating SO coupled BECs is time- 
dependent in the rotating frame, which is a considerably 
more complicated problem than the usual rotating BECs. 
Nevertheless, below we only consider the situation of the 
isotropic Rashba SO coupling, such that it in principle 
can be implemented as a time-independent problem in 
the rotating frame. 

The effect of rotation should be described by the stan- 
dard minimal substitution method as presented in Ref. 
[74j . The non- interacting part of the Hamiltonian is 

H = f d 3 r^(r) [^HftV + MXz x a — A) 2 — fi 



1 



+ V ext (r) - -hml{x 2 + y 2 )l XrO, 



(35) 



where the last term is the centrifugal potential due to 
rotation. 

Note that due to the presence of SO coupling, we 
should carefully distinguish the difference between me- 
chanical and canonical angular momenta. The mechani- 
cal one should be defined according to the minimal sub- 
stitution as 



ch = L z + M\(x<T x +y<T y ), 



(36) 



where L z is the usual canonical angular momentum. Ex- 
panding Eq. [35l it is equivalent to Eq. [T]plus the term 
of angular velocity Vt z coupling to L mech as 



d 3 n/^(r) 



L 



rnech 



Vv(r)- 



(37) 



Thus in the rotating frame, the effect of fl z is not only 
just fl z L z as usual, but also an extra effective radial Zee- 
man term as 



B R (r) = n z MXf. 



(38) 



Such a term is often missed in literatures. As will be 
shown below, it affects the ground state vortex configu- 
rations significantly and thus should not be overlooked. 

To make the model more adjustable, an external spa- 
tially dependent Zeeman field B ex — Br is intentionally 
introduced as 

Hb = —B [ d\ ^t(r) (xa x +ya y )^ </v(r), (39) 



which shares the same form as Eq. 1381 Experimentally, 
such a Zeeman field can be generated through coupling 
two spin components using two standing waves in the 
x and y-directions with a phase difference of ^. The 
corresponding Rabi coupling is 



n 



'{ sin(fc L x) + ism(k L y)^l(r)ij.^(f) + h.c. (40) 
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FIG. 6: The spin density distribution of the condensate of 
Fig. [5] (d). The projection of (a) in the xy-plane is shown 
as black vectors. A color map is used to illustrate the (cr z ) 
component. The unit of length for the figure is It- From Ref. 



In the region of \y\ <§; 2-x/kL, it reduces to the desired 
form of Eq. [35] with B — Vl'ki,. Such a term compensates 
the non-canonical part of the mechanical momentum in 
H rot , which renders the model adjustability in a wider 
range. 

B. SO coupled bosons in rotating traps 

Now we turn on interactions and obtain the ground 
state condensate numerically by solving the SO coupled 
GP equations which have been reduced into the dimen- 
sionless form as 

Ufa = T^+p(\^\ 2 + \^\ 2 )^, 

= f^ v + mil 2 + l^tl 2 )^, (41) 

where (ip^, tpi) are normalized according to the condition 
Jdr 2 ^ 2 + |^| 2 ) = 1. f is defined as 



T = -^l 2 T(dl + d 2 y ) + al T {-id y o x +id x a y ) 



^t(x 2 + y 2 ) - j-{-ixd y + iyd x ) 

a—(xa x + ycTy), 
It 



(42) 



where p = f2 z /u;; k = 7 + p; and 7 = B /(MuiX) is defined 
for the extra radial Zeeman field in Eq. 131)1 



1. The skyrmion lattice structure in the weak SO coupling 

Rich structures of vortex lattices appear in the case of 
weak SO coupling. In Fig. [5J the SO parameter is taken 
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FIG. 7: From left to right: the density and phase profiles of 
tp-f(r) and i>i(r) with parameter values of a = 4, j3 — 20, c = 
1, and p — 0.1. From (a)-(ft), 7 is taken as 0.5, 0.3, 0.1, —0.05, 
—0.25, —0.35, —0.6, and —0.7, respectively. The black arrow 
in each domain represents the local wavevector direction of 
the corresponding plane-wave state, which shows a clockwise 
or counter-clockwise configuration depending on the sign of 
7. For sufficiently large values of | , condensates distribute 
around a ring in space forming a giant vortex. The color 
scales for the density and phase distributions are the same as 
that in Fig. [5] The black circle with an arrow indicates the 
the direction of the circulation around the vortex core. The 
unit of length for the figures is It- From Ref. ff^ . 



as a = 0.5. Both the density and phase patterns for ipf(r) 
and 4>li r ) are depicted. Let us first consider the case of 
pure rotation of B = as shown in Fig. [5](d), i.e., 7 = 0. 
For ip^(r), its density distribution exhibits several dis- 
connected peaks. By contrast, the usual vortex lattices 
show disconnected low density vortex cores. The phase 
distribution V , t( r ) exhibits singular points around which 
phases wind 2tt. These singular points are squeezed out 
to the low density region near the edge. For ipi{r), its 
vortex cores are pinned by peaks of the density of ip^{r). 
Combine the distributions of Vf ( r ) an d 4>i(r) together, 
the condensate exhibits a skyrmion lattice configuration 
with the spin distribution (S) shown in Fig. [5] 

Turning on the external Zeeman field B ex of Eq. 1551 
changes the lattice configuration. If B ex || Br, the pa- 
rameter 7 > 0, and otherwise 7 < 0. For both 7 > 
and 7 < 0, if I7I is small, the skyrmion lattice struc- 
tures remain as depicted in Fig. 0(6, c, e, /). Increasing 
I7I further, the condensates of both spin components are 
pushed outwards, and distribute around a ring with a 
giant vortex core, as shown in Fig. 0(a) and (g). This 
ring is the location of potential minima shifted from the 
trap center by the Hb term. In all cases in Fig. \5\(a-g), 
the difference of vortex numbers between the spin- up and 
down components is one. This is a characteristic feature 
brought by the Rashba SO coupling. 



2. The domain wall structure in the strong SO coupling 

As discussed in section UlIBl in the case of strong SO 
coupling, if interactions are also strong, condensates are 
nearly suppositions of SO coupled plane- wave states sub- 
ject to the trap boundary condition. We consider the ef- 
fect of rotation in this case. Limited by the numeric con- 
vergence, only a small rotation angular velocity is con- 
sidered. The term of the external Zeeman field B ex of 
Eq. [35J is also applied, which enriches the structures of 
the condensates as shown in Fig. [7] 

The characteristic feature is that the condensate 
around the trap center is broken into several domains. 
Inside each domain, the condensate is approximated as 
a plane-wave state. Wavevectors are arranged such that 
the local spin polarizations are parallel to B ex in order 
to minimize the Zeeman energy. At small values of I7I as 
shown in Fig. [7](c, and d), a line of vortices appear at the 
boundary to separate two adjacent domains. Depending 
on the direction of B ex , the wavevectors inside domains 
can be clockwise or counterclockwise. For instance, when 
7 > 0, a clockwise configuration of these local wavevec- 
tors is selected in the condensate. As further increasing 
I7I, different domains connect together to form a giant 
vortex as shown in Fig. [7] (a, g, h). Both spin compo- 
nents overlap with each other, and distribute around a 
ring with the radius of a|7|iy. The spin textures lie along 
the radial direction to minimize the magnetic energy. 
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V. MAGNETIC PHASES OF SO COUPLED 
BOSONS IN OPTICAL LATTICES 

It is natural to further consider the SO coupling effect 
in optical lattices, in particular in the Mott-insulating 
states^— . Due to SO coupling, hopping amplitudes 
are spin-dependent whose values vary non-monotonically 
as increasing SO coupling strength. The spin-dependent 
hopping leads to the Dzyaloshinsky-Moriya (DM) type 
spin exchange models in the Mott-insulating states, 
which results in rich spin ordering patterns. 



Tight-binding approximation and band 
parameters 



(a) 

(Er) 
0.03 

0.02 

0.01 

0.00 

-0.01 - 



V 



\ 



\ 
\ 



X /' 



0.4 0.6 0.8 

KJK 

SO 



FIG. 8: The dependence of the spin-independent hopping 
integral t and the spin-dependent one A v.s. the SO coupling 
strength k ao /ko. The optical potential depth is Vo = 8E r . 
From Ref. U^. 

The tight-binding model in the square optical lattice 
will be derived below for two-component bosons with syn- 
thetic SO coupling. The single-particle Hamiltonian in 
the continuum is defined as in Eq. [T] by replacing the 
trapping potential with the periodic lattice potential as 



V(x, y) = —Vq cos 2 kox + cos 2 k y 



(43) 



where fco = 2-7r/Ao and the lattice constant a = Ao/2. 
The recoil energy is defined as E r — h 2 k^/ (2M). For 
later convenience, the relative strength of SO coupling 
is quantified by the dimensionless parameter k so /ko with 
k so = MX R . 

A tight binding Hamiltonian for the lowest orbital 
band with SO coupling can be written as^ 



H = — 



^ i J; cr " 

{ij) ,c 



(44) 



where only the nearest neighbor hoppings are included. 
Uj-o-o-' can be decomposed into spin-independent and 



spin-dependent components based on the following sym- 
metry analysis 



t 



= t-\-i\i 



(45) 



The coefficient of spin-dependent hopping is purely imag- 
inary as a requirement of the TR symmetry. The optical 
lattice with the Rashba SO coupling possesses the reflec- 
tion symmetry with respect to the vertical plane passing 
bonds along the x and y-directions. These reflection sym- 
metries require that Xi,i+e m || e y , and Xi t i+e y \\ —e x . The 
four-fold rotation symmetry requires that A^+e^ ■ &y = 
— ^i,i+e H • &x- All these symmetry properties together 
constraint the spin-dependent hopping up to a single pa- 
rameter A as 



A, 



= Ae 



Xi,i-\-e y 



(46) 



The band structure parameters t and A are related to 
the overlap integrals of the onsite SO coupled Wannier 
functions in neighboring sites. In the case of a deep lat- 
tice, each site can be approximated by a local harmonic 
potential. The lowest energy Wannier states ipj z= ±i 
are a pair of Kramer doublets as presented in Eq. [33] 
The radial wavefunctions of the Wannier states f(r) and 
g(r) exhibit the Fricdcl-typc oscillations as explained in 
Sec. IIII Al Thus naturally t and A should also exhibit 
such oscillations as increasing the SO coupling parame- 
ter k so /ko. This feature is numerically confirmed using 
the following method. The tight-binding band spectra 
can be calculated easily as 



E±(k) = e(k) ± 2Ai/sin 2 k x + sin 2 k 



1 y 



(47) 



where e(k) = —2t(cosk x + cosk y ). On the other hand, 
the band spectra can be calculated directly from the con- 
tinuum model with the lattice potential Eq. [43] by using 
the basis of plane- waves. By fitting these spectra using 
Eq. 27] the values of t and A are obtained and are plotted 
in Fig. [5] Both t and A oscillate as increasing fc so /fc , and 
the amplitudes of their overall envelops decay. For the 
spectra of Eq. 27] the square lattice breaks the rotational 
symmetry down into the 4-fold one, thus the degeneracy 
of the Rashba ring is lifted. The lower band has the 4-fold 
degenerate minima located at Q = (±fc, ±fc) with 



tan 



1 A 



(48) 



It should be pointed out that these lowest local Wan- 
nier states are the eigenstates of the on-site total angular 
momentum with j z — ±i. Therefore, the eigen-bases 
defined by (bif,bii) T should be those of j z not a z . In 
the case of strong SO coupling k so > fco, the angular 
momenta of these Wannier states nearly come from the 
orbital angular momentum moment, while their spin mo- 
ments are nearly average to zero. 

If in the limit of very strong SO coupling such that 
k so 3> fco, Landau level quantization effects appear 
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within each site. Many states with different values of 
j z are nearly degenerate as presented in Sec. Ill Bl In 
this case, a single band model Eq. |44] fails even in the 
case of the deep lattice. It is justified only in the case 
that k so /k < 1 in which the lowest Wannier states are 
separated from others. 

B. Magnetic properties in the Mott-insulating 
state 

We consider the spin physics in the Mott-insulating 
phase of Eq. 02] In the simplest case, there is one par- 
ticle per site with a two-fold degenerate Kramer dou- 
blet. The low energy superexchange Hamiltonian can 
be constructed using the Schrieffer-Wolf transformation, 
which shows the Dzyaloshinsky-Moriya (DM) type ex- 
change due to SO couplin g 103 ! 104 as 

i 

and 

Hi^i+e^ = —J\Si ■ >Sj+e tl — J^di.i+e^ 1 (Si X Sj+g^ ) 

+ J2[Si ■ S l+ g tl — 2(Si ■ di t i+g tl )(Si+£ fi ■ d%>i+&p)u 

(50) 

where e^(fi — x,y) are the unit vectors along the x and 
y-directions, respectively; J\ = 4t 2 /U, J\2 = At\/U, and 
J 2 = A\ 2 /U. The DM vectors are defined as df j+g,. — &y 
and dij+s = —e x which are perpendicular to each other. 
This is similar to the case of the high T c cuprate su- 
perconductors such as YBa2Cu30 u 105 ' 106 . Consequently, 
these DM vectors in Eq. [50] cannot be removed by gauge 
transformations, or, equivalently by varying local spin 
axes. This brings frustrations to magnetic properties. 
To obtain a qualitative understanding, two different lim- 
its of |A| <C \t\ and |A| 3> \t\ will be considered. 

In the absence of SO coupling, i.e., A = 0, the system 
is in the ferromagnetic state. If A is small, the J2-term 
brings the easy plane anisotropy which prefers spin mo- 
ments lie in the xy-plane. The DM-vector further in- 
duces spin spiraling at a finite wavevector, which can be 
shown by calculating the spin-wave spectra around the 
variational ground states that spin moments lie along the 
high symmetry line of diagonal directions, say, [110]. The 
Holstein-Primakoff transformation is employed to trans- 
form Eq. (|49p into the magnon Hamiltonian, 

T f/ <->/! sin 26* t 
H mg = - Jo ((cos 26» - % — ) a ! a i+ex 

+(cos29 + i h -^-)a\a. l+ey +h.cj, (51) 

where a 1 is the creation operator for magnons deviat- 
ing from the [110]-direction; 8 = arctan(A/t) as defined 
above. We only keep quadric terms and ignore the terms 



proportional to sin 2 6 since A/t < 1. In momentum 
space, its spectra can be diagonalized as 

e(k) — — 2 Jo| cos 2(9(cos k x + cos k y ) 

+ —= sin 26 , (sin fca; — sinfc y )|, (52) 

whose minima are located at Qm = (2fc, —2k) with the 
value of k given in Eq. 05] This indicates that the ground 
state exhibits a spin spiral order along the direction per- 
pendicular to the quantized axis in the spin-wave analy- 
sis. 

Interestingly, in the opposite limit of |A/f| 3> 1, Eq. 
[50] can be related to that of |A/i| <C 1 through a duality 
transformation. On site i with the coordinates (i Xl i y ) 1 

Si is transformed into 

su -+ (-if^u- (53) 

Si still maintains the spin commutation relation. Under 
this transformation, the Ji-term transforms into the J 2- 
term and vice versa, and the Ji2-term is invariant. Thus 
this dual transformation indicates that there is a one- 
to-one correspondence between the ^-dominant phase 
(|A/i| ^> 1) and that of Ji with \X/t\ <C 1 which has been 
analyzed above. 

In the regime of intermediate values A/t, a rich phase 
diagram with different spin patterns appears. Clas- 
sical Monte Carlo simulations have been employed to 
calculate the ground state phase diagram in current 
literatures^—. Various patterns have been found as 
a result of competition among ferromagnetic exchange, 
easy-plane anisotropy, and the DM effect induced spin 
spirals. These include the ferromagnetic, antiferromag- 
netic, spiral, stripes, and vortex crystal orderings. Fur- 
thermore, the superfluid-insulator transition for SO cou- 
pled bosons has also been studied in Ref. [79| . 

VI. CONCLUSIONS 

We have reviewed unconventional BECs with SO 
coupling whose condensate wavefunctions are complex- 
valued and are thus beyond the framework of the "no- 
node" theorem. Even at the single particle level, 
the spectra in harmonic traps exhibit the structure of 
Landau-level like quantization induced by SO couplings. 
Their energy dispersion is nearly flat with respect to an- 
gular momentum in the case of strong SO coupling, and 
exhibit the Z2-type topology. The interacting conden- 
sates exhibit topologically non-trivial configurations. In 
the 2D Rashba case, the spin density distributions are 
characterized by the skyrmion type textures. The 3D 
Weyl SO coupling induces the topological phase defects in 
the quaternionic phase space, and the corresponding spin 
density distributions are also non-trivial carrying non- 
zero values of the Hopf- invariant. In rotating traps, the 
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condensate configurations are changed by vorticity which 
results in a variety of structures including skyrmion lat- 
tices, giant vortices, multi-domains of plane-waves. In 
the strongly correlated Mott-insulating states, SO cou- 
pling exhibits in the DM exchange interactions in the 
quantum magnetism. The research of the novel states of 
SO coupled bosons is still in the early stage. In partic- 
ular, the effect of SO couplings in the strong correlation 
regime is still a largely unexplored held. We expect that 
further exciting progress on the novel states of SO cou- 
pled bosons will appear in the near future. 
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